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We address the joint estimation of the two defining parameters of a displacement operation in 
phase space. In a measurement scheme based on a Gaussian probe field and two homodyne detectors, 
it is shown that both conjugated parameters can be measured below the standard quantum limit 
when the probe field is entangled. We derive the most informative Cramer-Rao bound, providing 
the theoretical benchmark on the estimation and observe that our scheme is nearly optimal for a 
wide parameter range characterizing the probe field. We discuss the role of the entanglement as 
well as the relation between our measurement strategy and the generalized uncertainty relations. 

PACS numbers: 42.50.Xa, 03.65.Ta 



I. INTRODUCTION 

One of the most promising avenues for quantum tech- 
nology is the use of quantum resources to improve the 
sensitivity in the estimation of (not directly observable) 
relevant physical parameters, e.g. for applications in 
metrology and sensing [TJ |2]. While the estimation of 
a single parameter has been extensively studied from 
both theoretical and experimental viewpoints p] [3] , the 
joint estimation of multiple parameters has not received 
enough attention (notable contributions are [3HL3]). Here 
we take a further step forward in this direction by pro- 
viding optimal quantum-enhanced strategies to estimate 
the two conjugated parameters characterizing a paradig- 
matic and ubiquitous quantum operation, phase-space 
displacement. Quantum states of single-mode bosonic 
CV systems can be described by quasiprobability distri- 
butions on a two-dimensional real phase space The 
operation of displacing a state by a phase space vector 
(lo,Po) is represented by the Weyl displacement operator 



D(qo,po) = exp(ip <7 - iqop) , 



(1) 



where q and p are the two quadrature operators satisfying 
the canonical commutation relation [q,p] — it. A rele- 
vant question is, starting with a reference state go that 
undergoes an unknown displacement, how accurately can 
we jointly estimate the two conjugate parameters go an d 
Po of the displacement operator with a measurement on 
the displaced state g = D(qo,po)goD 1i (qo,Po)^ One pos- 
sibility is certainly to use coherent states as initial probe 
states, followed by heterodyne detection as a measure- 
ment strategy at the output of the displacement transfor- 
mation 0] . This has in fact been the standard technique 
to estimate and measure displacement, as it naturally 
complies with the generalized Heisenberg uncertainty re- 
lation [17] . On the other hand, one may ask whether 
entanglement, in the form of Einstein-Podolsky-Rosen 
(EPR) correlations [T8|, could lead to a better estima- 
tion precision for this task, as it was suggested in |19j . 



as it happens in the case of phase estimation in quantum 
metrology [T] and for the joint estimation of the tem- 
perature and coupling constant in a bosonic dissipative 
channel |14) . 

In this paper we propose and analyze a measurement 
scheme where CV entanglement is used to improve the 
estimation precision for this particular and relevant prob- 
lem. For this purpose, we limit the analysis to Gaus- 
sian states and operations [30] . An introduction on LQE 
for multiple parameters will be given in Sec. [TTJ while 
in Sec. [Hi] we will derive the bounds for the displace- 
ment estimation by considering single- and two-mode 
Gaussian states. In Sec. |IV| we show a simple mea- 
surement scheme involving two-mode squeezed thermal 
states which achieves the ultimate bound for different 
values of squeezing and thermal photons and beats the 
standard quantum limit for this kind of estimation. We 
also discuss how the performances and the bounds change 
when one has a-priori information about the parameters 
to be estimated, and when the displacement operation 
presentes an inner uncertainty. In Sec. |V]we discuss the 
role of entanglement, showing how it is always necessary 
to beat the classical optimal strategy, and also necessary 



for symmetric probe states. In Sec. VI we discuss the re- 



lationship between this multi-parameter estimation and 
the generalized uncertainty relations, and finally we end 
the paper with some concluding remarks. 



II. MULTI-PARAMETER LOCAL ESTIMATION 
THEORY 

How can we know that a given measurement scheme, 
able to estimate certain parameters, is optimal? Is it 
possible to estimate these parameters with a better pre- 
cision? In order to answer this question, we make use 
of tools derived from local quantum estimation (LQE) 
theory [2l l2lT - f23] . The purpose of LQE theory is indeed 
to determine the ultimate precision achievable and the 
corresponding optimal measurement for the estimation 
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of parameters characterizing a physical quantum system. 
In particular it has been applied for the estimation of dif- 
ferent quantities, including the quantum phase of a har- 
monic oscillator [2H [55] , CV Gaussian unitary parame- 
ters |26) , the amount of quantum correlations of bipartite 
quantum states j27j[28], and the coupling constants of dif- 
ferent kinds of interactions [T31 RHI I2M3T] . While most 
studies so far have been limited to single parameter esti- 
mation, the case of multiple parameters is more complex 
as different bounds can be derived for the same setting. 
Moreover, as we will point out later, these bounds are 
not always achievable, in particular when one deals with 
conjugate variables for which a Heisenberg-type uncer- 
tainty relation applies. 

Let us start by considering the general case, that is a 
family of quantum states g z which depend on a set of d 
different parameters z = {z^}, fi = 1, ••• ,d. One can 
define the so called Symmetric Logarithmic Derivative 
(SLD) and Right Logarithmic Derivative (RLD) opera- 
tors for each one of the parameters involved, respectively 
as 



dg 
dg 



Lf ] g + gL^ ] 



SLD. 
RLD. 



Then one can define the two matrices 



t(S)t(S) 



H™ = tr 



V = tr[ & 4 fl) 4 fl)t ] 



T {S) T (S) 



(2) 
(3) 

(4) 
(5) 



By defining the covariance matrix elements V^z)^ ~ 
E[z^z v ] — E[z^]E[z v ) and a weight (positive definite) ma- 
trix G, two different Cramer-Rao bounds holds 



RLD bound Br on the sum of the variances are in gen- 
eral achievable 5J. The first one could not be achiev- 
able because it corresponds to the bound obtained by 
measuring optimally and simultaneously each single pa- 
rameter and this is not possible when the optimal mea- 
surements do not commute. At the same time the RLD 
bound could not be achievable because the optimal esti- 
mator does not always correspond to a proper quantum 
measurement (that is a proper positive operator valued 
measure). Moreover which one of these bounds is more 
informative, that is which one is higher and then tighter, 
depends strongly on the estimation problem considered. 
One can then define the most informative Cramer-Rao 
bound, 

B M i = max{B s , B R } , 
obtaining the single inequality 

Bmi 



E Var (^)> 



M 



A. Cramer-Rao bounds with a-priori information 

Similar bounds can be obtained in the case where one 
has a certain a-priori information regarding the distri- 
bution on the parameters one wants to estimate. Let 
us assume that the a-priori information is described by 
a probability distribution V pr \ or {z). One can define a 
Fisher-information matrix of the a-priori distribution as 



A^„ = / dzV prlor {z) 



9l0g7 , prior(z)\ (d\0gV 

prior ) 



dz u 



dz v 



(10) 



tr[GV] > itrlGtH)" 1 ] 



tr[GV] > 



trfGRcfJ- 1 )] + tr[|GIm(J~ 1 ) 
M 



(6) 
(7) 



where trL4] is the trace operation on a finite dimensional 
matrix A and M is the number of measurements per- 
formed. We observe that if we choose G = 1, we obtain 
the two bounds on the sum of the variances of the pa- 
rameters involved 



1 

M 



trfH- 1 ] 



trfRc^- 1 )] +tr[|Im(J" 1 ) 
M 



(8) 
(9) 



The matrices H and J arc called respectively the Sym- 
metric Logarithmic Derivative (SLD) [22] and Right Log- 
arithmic Derivative (RLD) [U [BHH] quantum Fisher in- 
formation matrices. Neither the SLD bound Bs, nor the 



The new Cramer-Rao bounds that take into account this 
a-priori information, will read 



trlGVl^^tr^H + A)- 1 ], 

tr[GV] > — (tr[GRc((J + A)" 1 )] + tr[|GIm((J + A^ 1 

1V1 



(11) 

A)- 1 
(12) 



and, for G = 1, 

£Var(z,)> B ^ : = ^tr[(H + Ar l ] 
£Vax(, M )>^ 



^-(trlReaj + A)- 1 )] 
-trOMCJ + A)- 1 )!]) . 



(13) 



(14) 



Here A denotes a vector of parameters characterizing the 
prior information at our disposal. 
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B. Evaluation of the RLD Fisher information 

In the following we give some details about the deriva- 
tion of the RLD Fisher information when the RLD op- 
erator cannot be evaluated directly. Let us suppose that 
the derivative with respect to every parameter has the 
following form 



dg 



BuQ- 



O) 



(15) 



To obtain the RLD operator L M = 
in Eq. ^ we have to find the operator L M 
(assuming that g _1 exists) 

B M g = Ql$ . 



as defined 



(b) 



such that 

(16) 
(17) 



the param eter s /i = {qo,Po}, we obtain formulas resem- 

-iG/j,. Then, starting 



bling Eq. (15 



where L 



(a) 
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from Eq. (T8), and by observing that 



(22) 
(23) 



we can express, after some algebra, the elements of the 
Fisher information matrix in terms of the generators of 
the displacement as 

J^u = tr[G v $GnQo : ] + tr[£> G„G M ] - 2tr[g G M G„] . 

We notice that the Fisher matrices do not depend on the 
values of the parameters to be estimated and that the 
only elements that are involved are the probe state and 
the generators of the two transformations. 



Then, after some algebra, we can express the elements of 
the RLD Fisher information matrix as 



= tr[gLW(LW)t] 
+ tr[ e (L^B v ] 



-tv[gBlL^}+ 
MBufBlQ- 1 ] 



(18) 



III. CRAMER-RAO BOUNDS FOR 
DISPLACEMENT ESTIMATION 

Let us now consider a generic probe state (pure or 
mixed) go, which is displaced by the operator D(qo,po) 

to the state g = D(qo,po)goDHlo,Po)- In the following 
we derive explicit formulas for the SLD and RLD Fisher 
information matrices. Let us start by considering the 
SLD Fisher information for a given probe state whose 
diagonal form reads go — ^2 n Pn\4>n){4>n\- One proves 
that in our case the SLD operator in Eq. pi) satisfies the 

property = D(q , p )C f _ l D^ (q 0l po), where 



£„ = 2i (Gu) nm — Pm \<t> n ){4>-n 

-V Pn + Pm 



(19) 



with = {<?o,Po}, (G^nm = (</>n|G M |</> m ) and where 
G qo = p, G po — —q are the generators of the two orthog- 
onal displacements. Then the SLD Fisher information 
matrix elements read 



H 



(20) 



25> 



Ps - Pt 
Ps +Pt 



((G M ) s t(G 1/ )f S + (Gjy) stiG/jjts) 

(21) 



Let us consider now more in detail the case of the RLD 
Fisher information. By differentiating g with respect to 



A. Most informative bounds for single- and 
two-mode probe states 

The most general single-mode Gaussian state with zero 
initial displacement can be written as go — S(r)vNS(ry 
where = jf+i^2 n j^i\n}(n\ is a thermal state and 
S(r) — exp{— \[a^ 2 — a 2 )} is the single-mode squeez- 
ing operator. Notice that every single-mode squeezed 
state evolving in a noisy dissipative channel can always 
be written in this form, which makes this treatment im- 
portant for actual implementations [38 . In this case 
the two bounds Bg and Br are evaluated and the most- 
informative for the single-mode case b[J| is found to be 
equal to the RLD bound, yielding 



B Mi( r > N ) = (27V + 1) cosh 2r + l 



(24) 



For zero squeezing the results obtained by Yuen and Lax 
is recovered |3] . Moreover one can verify that this bound 
is achieved for any value of squeezing and thermal pho- 
tons by performing an heterodyne measurement. In gen- 
eral we observe that the bound grows with N and r. It 
is thus clear that single-mode squeezing is not useful for 
displacement estimation and the optimal measurement 
setup involving single-mode Gaussian probe states and 
heterodyne detection corresponds to using the vacuum 
(or any coherent state) as a probe field. The correspond- 
ing bound is denoted by 



B 



sql 



B«(0,0)=2 



(25) 



as to the standard quantum limit (SQL). We note that 
the SQL does not depend at all on the mean energy of 
the probe coherent state: by increasing the mean photon 
number of the coherent states one does not obtain any 
enhancement in the estimation precision. 
Let us focus now on the more interesting two-mode case, 
where the displacement operator is applied only on one 
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FIG. 1. (Color on-line) Measurement scheme for the esti- 
mation of the displacement given a two-mode squeezed probe 
state. After the displacement operation, the modes are mixed 
in a balanced beam splitter and then orthogonal homodyne 
measurements are performed on the output modes. 



of the two-modes. The probe state corresponds to a two- 
mode squeezed thermal state, which is an archetype of 
the (possibly noisy) Gaussian entangled states: 



Qo 



S 2 (r){v N ®VN)SUr) 



(26) 



where S2O") — exp{r(aW — ab)} is the two-mode squeez- 
ing operator. The two bounds can also be straightfor- 
wardly evaluated, obtaining 



B^(r,N) 
B^(r,N) 



2N + 1 
cosh 2r 

AN(1 + N) 
(2iV+l)cosh2r- I' 



(27) 
(28) 



Both are increasing functions of the average number 
of thermal photons N and decreasing functions of the 
squeezing parameter r (and thus of the entanglement of 
the probe state). In this case which bound is the most 
informative depends on the actual values of r and N. 
Comparing Eqs (36) and (37), when cosh(2r) < 2N + 1, 
Bs < Br. Thus we define a threshold value for the 
squeezing as r ths = ~ cosh~ 1 (2A r + 1), and the most in- 
formative bound reads 



Bf\r,N) 



for r < r t hs 
for r > r ths . 



(29) 



We notice that for N = the most informative bound 
coincides with the SLD bound, while if we increase the 
value of N and for small values of the squeezing param- 
eter r, the most informative bound turns out to be the 

RLD bound. By inspecting the most informative bound 

(2) 

°mv we n °tice that for different values of the parame- 
ters the bound is smaller than the SQL bound B sq i. One 
may then wonder if by using entangled probe states one 
can achieve a better result; in the next section we indeed 
present a simple measurement scheme, outperforming the 
classical single-mode strategy. 



IV. NEAR-OPTIMAL MEASUREMENT 
SCHEME 

As pointed out in the previous section, if we consider 
coherent states as probe states and then, after the dis- 
placement operation, we perform a heterodyne measure- 
ment, we achieve the SQL B sq i. On the other hand, the 
bounds obtained for entangled probe states suggest that 
the SQL can be in principle overcome. Indeed we now il- 
lustrate the two-mode measurement scheme able to beat 
bound B sq i and to achieve the optimality for different val- 
ues of the parameters characterizing the probe state. The 
scheme is pictured in Fig. [I] It clearly resembles the CV 
version of the dense coding protocol [39] and was already 
suggested for the estimation of displacement [19 . The 
probe state corresponds to a two-mode squeezed thermal 
state ( 26 ) . The displacement operator is applied on one 



mode after which the two modes are mixed at a balanced 
beam splitter. Then the output fields of the beam splitter 
are described by the density operator 



d = U hs gU hs = g\ <g> g' 2 



(30) 



where Ubs = exp{^(af>^ — a^&)} is the beam splitter oper- 
ator [TB], g = D(q ,p )g & (qo,Po) is the state after the 
displacement, and 

Q'i = D(qip' )S(r)iy N S\r)D^q^p f) (31) 
g' 2 = D(q' ,p' )S(-r)v N SH-r)D\ q ' ,p' ) (32) 

with q' = ^ , p' = ^ [ID] . The output state is a tensor 
product of two states squeezed in orthogonal directions 
and both displaced by the rescaled values q' and p' a . One 
performs a homodyne measurement of the quadrature p 
on the state g[ and of the quadrature q on g' 2 , obtaining 
respectively the parameter values go an d Po- As the states 
are squeezed in orthogonal directions, the two variances 
approach exponentially to zero by increasing the squeez- 
ing parameter r as Var(g ) = Var(p ) = (2A^ + l)e~ 2r : 
the higher is the squeezing the more precise is the esti- 
mation. The sum of the two variances is 

Var( go ) + Var(p ) = 2 (2A^ + 1) e~ 2r > B^j(r, N) (33) 

One can observe that we obtain for the two-parameter 
estimation the same optimal scaling in terms of the de- 
gree of squeezing, as the one obtained for the single- 
parameter displacement estimation in [41] ; in particu- 
lar for a pure two-mode squeezed state (N — 0), one 
achieves for large squeezing the Heisenberg limit scaling 
l/N, where N = sinh 2 r denotes the mean number of 



photons. Comparing Eq. (33) with Eq. (25), it is clear 



that this scheme can outperform the single-mode strat- 
egy. For TV = 0, as long as squeezing is non-zero, we can 
estimate the parameters better than the SQL suggests. 
For N 0, if the field exhibits two-mode squeezing, that 
is, if it is squeezed stronger than the following threshold 



r > r sql (N) 



1 



In 1 1 + 4N 



AN 2 ) 



(34) 



5 




FIG. 2. (color on-line) Renormalized difference D(r, N) 
between the sum of the variances for the estimation of dis- 
placement with the double-homodyne scheme and the most 
informative bound B^j] (r, N) , as a function of the squeez- 
ing parameter r and for different values of thermal photons: 
continuous-red line N = 0; dashed-green line N = 0.5; dotted- 
blue line N = 2. 



where G^,a 2 {%) denotes a Gaussian distribution, centered 
at [i and with variance a 2 . The objective is to minimize 
the average precision one gets on the estimation of these 
random parameters; we can then evaluate the bounds in 
Eqs. (13) and (14), for different input states and strate- 
gies. By considering coherent states as the input, the 
most-informative RLD bound will be equal to 



Bsql(A) := 



2A 2 



1 + A 2 ' 



(35) 



If we rather consider a two-mode squeezed thermal state 
f?o = S2(r)i'N <8> vnS\(t) as the probe, we obtain the 
following bounds: 



B g 2) (r,7V,A) 
B$\r,N,A) 



2(27V + 1)A 2 



(36) 



2N + 1 + 2A 2 cosh 2r 

AN (I + N)A 2 
2N(1 + N) + A 2 [(27V + 1) cosh2r - 1] ' 

(37) 



we can beat the SQL. 

If we now compare the obtained results to the most- 
informative bounds derived in the previous section, we 
observe that the sum of the variances E(r, N) = Var(g ) + 
Var(po) in Eq. (33) is, as expected, bounded from below 



(2) 

by Bj^j(r, N). However one may wonder if, in some range 
of parameters, the scheme becomes (nearly) optimal, that 
is if the bound is almost saturated. In Fig. [2]we plot the 
quantity 



D(r,N) 



E(r,N)-B™{r,N) 
<(r,N) 



for different values of N and as a function of the squeezing 
r. One observes that by increasing the squeezing param- 
eter r, our scheme is optimal with D(r, N) ~ 0. For lower 
values of r, we notice that D(r, N) is not always mono- 
tonically decreasing; this is because the most informative 
bound changes between the RLD bound (r, TV) and 

the SLD bound B g '(r,N), as explained before. As re- 
marked before, the measurement scheme we present re- 
sembles the CV version of the dense coding protocol [35] ; 
however, while the dense coding protocol requires more 
than 4 dB of squeezing to outperform single-mode strate- 
gies, our estimation strategy outperforms the SQL for 
any value of two-mode squeezing at the input. 



A. Estimation with a-priori information 

Let us consider now the case where we have some 
prior distribution on the parameters we want to estimate. 
In particular, for the sake of simplicity we consider the 
parameters taken randomly with the following a-priori 
probability distribution: 

^prior(go,Po) = £o,A(go) 0O,a(Po), 



All these bounds decrease with the decreasing A. We 
note that all the bounds discussed before can be re- 
obtained by taking the limit of flat a-priori distribution 
(A — > oo ). If one fixes the value of A, one can de- 



fine the most-informative bound B^j = m&x{B { g' , B^" 1 }, 
and compare it with the corresponding SQL bound. One 
could then ask which is the optimal measurement strat- 
egy, and if the precision obtained saturates the most in- 
formative bounds for different possible probe states. Let 
us start by considering input coherent states: as proved 
in [3], the optimal cheating strategy corresponds to mul- 
tiply the heterodyne outcomes by a factor 



(2) R (2) 



K c = 



A 2 



1 



(38) 



It is indeed easy to check that with this choice the ob- 
tained averaged variances are equal to the SQL limit 
Bsql(A) derived in Eq. ( [35] ) (notice that this is also 
the optimal choice used in |421 14"3"] to derive the classi- 
cal benchmark for teleportation of coherent states). Let 
us consider the general case where, for given values of 
qo and po, the variances obtained with a certain mea- 
surement strategy are equal and do not depend on the 
parameters themselves, that is 

Var(g ) = Var(po) := Var . 

One can prove that the scaling factor that minimizes the 
average sum of the variances is equal to 



A 2 

K ■ = 

1 1 mm -.j . ao' 

Var + A z 



(39) 



and the obtained result is 

( Var^(go) + Var^(po) ) 



2Var A 2 
Var + A 2 
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where (•) denotes the average on the a-priori distribution. 
This is also the case for the two-mode squeezed thermal 
states considered before. Of course the scaling factor in 
this case depends on the probe state parameters, since 



Var = (2N + l)e 



-2r 



(40) 



If this information is not available, one can always adopt 
the coherent states optimal strategy and use the scaling 
factor K c in Eq. (38), which does not depend on the 



input state, obtaining 

A 2 f1 4- A 2 Var \ 
( Var^o) + Var* (p ) ) - (1 + A 2 )2 0> ■ (41) 

The different results for two-mode squeezed thermal 




0.5 1.0 1.5 2.0 2.5 3.0 



0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 3. (color online) Dashed-blue line: average sum of the 
variances for two-mode squeezed thermal probe states with 
N — 1 mean thermal photons and by adopting the optimal 
scaling factor Kmi n . Dotted-red line: average sum of the vari- 
ances for two-mode squeezed thermal probe states with N = 1 
mean thermal photons and by adopting the coherent state 
scaling factor K c . Solid-black line: most-informative bound 
for two mode squeezed thermal state (with N = 1). Solid- 
grey line: standard quantum limit Bsql- All the plots are 
functions of the squeezing parameter r and for different val- 
ues of the a-priori uncertainty A. Top-left: A = 1; top-right: 
A = 2; bottom-left: A = 3; bottom-right: A = 5. 

states are shown in Fig. [3] one observes that by using 
the scaling factor K m i n , the estimation strategy is nearly 
optimal, that is the most-informative bound is saturated 
for a wide range of parameters of the probe state, and for 
different values of the a-priori uncertainty A. One also 
observes that by using the simpler scaling factor K c , one 
still beats the SQL limit by increasing the squeezing pa- 
rameter; in particular, for zero thermal photons (N = 0), 
the entangled assisted strategy always beat the SQL for 
any value of the squeezing parameter r. On the other 
hand this strategy is far to be optimal for low values of 
A and for large values of the squeezing parameter r. 



B. Estimation of imperfect displacement 
operations 

Let us consider the case where the displacement opera- 
tion is imperfect. We thus have an additional uncertainty 
on the parameters we want to estimate. We assume that 
the two corresponding values are distributed according 
to a certain probability distribution V en (q / ,p') which has 
mean values qo and Pq. The output state, after the dis- 
placement operation, can thus be written as 

q= J dq'dp' V err (q',p')D(q',p')eDH<l',P')- (42) 

For the sake of simplicity, let us consider that the er- 
ror probability is a product of two Gaussian independent 
probability distributions, i.e. 

T , en(q',p') = G qa A 2 M')G P oA 2 p (p')- 

Using our entanglement-assited estimation strategy, we 
obtain the following result for the variances of the esti- 
mated paramaters: 

Var(go) + Var(p ) = 2(2N + 1) e~ 2r + A 2 , + A 2 ,. (43) 

It is clear that the additional uncertainties are simply 
added to the previous results, giving, as expected, a worse 
performance in terms of estimation precision. 



V. THE ROLE OF ENTANGLEMENT 

In our measurement scheme, we make use of entan- 
gled Gaussian states showing EPR correlations, as probe 
states. One may then ask whether the entanglement 
of these states is necessary (or even sufficient) to beat 
the SQL bound obtained by means of the single-mode 
strategy. For this purpose, we consider a generic two- 
mode Gaussian state, without local squeezing, i.e. with 
((Aqi) 2 ) — ((Api) 2 ), i = 1,2. This is a reasonable choice 
because we know that local squeezing does not help in 
our scheme. Given a generic two-mode quantum state, 
the corresponding quadrature operators cji and pi and an 
arbitrary (nonzero) real number a, if we define the oper- 
ators u and v as u — \a\q\ + ^(j2 , v = \a\p\ — , 
Duan et al. 



proved that, the condition 



((Aw) 2 ) + ((At;) 2 ) < a 2 + 



1 



(44) 



is a sufficient condition for inseparability. 
One can easily notice that, the inseparability condi- 
tion is the same as V&r(qo) + Var(po) < B sq i, assum- 
ing a = 1, which gives the lowest bound in Eq. (44 1 



This clearly shows that the entanglement of the probe 
state is a necessary condition if we are to beat the SQL 
obtained using coherent states and heterodyne measure- 
ments. Moreover, for symmetric states, such as the two- 
mode squeezed thermal state g = S2{r){vN ® v N)S\(r), 



7 



it is proved that the condition (44) with a = 1 is a nec 



essary and sufficient condition for inseparability |44j . As 
a consequence, for this class of states, entanglement is 
not only necessary but also sufficient to beat the SQL. 
It is also straightforward to find a counterexample in 
order to prove that in the asymmetric case, entangle- 
ment is only necessary but not sufficient. Let us con- 
sider an asymmetric two-mode squeezed thermal state 
Q = S 2 {r)v Nl ®VN 2 Sl{r), with N\ ^ N 2 ; if we set N x = 
the state is always entangled for r ^ 0, but to beat the 
SQL on the estimation of displacement, one can show 
that N2 has to be moderately low (one can derive the 
threshold value as a function of squeezing parameter r). 
It is worth stressing the fact that the state must be entan- 
gled before the application of the displacement operator. 
If we consider the case where a two-mode squeezer is ap- 
plied after the action of the displacement operator on a 
thermal state vn, no enhancement in the precision esti- 
mation can be achieved. Here, this squeezing operation 
can be thought as a part of the measurement process. 
The ultimate precision in this case coincides with the re- 
sults described for single-mode states and has to comply 
with the SQL. This result is in fact related to the secu- 
rity of the CV quantum key distribution protocol with 
coherent states |4"B"1. 



VI. PARAMETER ESTIMATION AND 
UNCERTAINTY RELATIONS 

We have observed that it is possible to measure the 
two conjugate parameters below the SQL. This seems 
to contradict with the (generalized) Heisenberg uncer- 
tainty relations. Nevertheless, if one looks carefully at 
the setup, one notices that the fundamental uncertainty 
relations are never violated: the variances corresponding 
to the true quantum quadrature operators qi and p$, on 
each mode involved and at every step of the measurement 
setup always satisfy the uncertainty relation, as it ought 
to be. The generalized uncertainty relations derived in 
|17j show that an inherent and unavoidable extra noise 
has to be taken into account if one wants to estimate two 
conjugate parameters by means of a joint measurement. 
However that analysis did not take into account the pos- 
sibility of having a two-mode entangled state as the initial 
probe as described in the previous scheme |46j . In fact, 
in our setup the pre-existent entanglement is exploited 
in order to perform precise measurements on different 
modes, and thus on commuting observables. Specifically, 
if we consider the product of the corresponding variances 
on the estimation of the parameters qo and po, we are 



led to conclude that the generalized uncertainty relation 
seems to be violated when 



Var(g )Var(p ) < 1. 



(45) 



If Var(go) = Var(p ), as it is always the case by consider- 
ing go = S2(r)vN(&i>NS\(r) as a probe state and our mea- 
surement setting, one can clearly observe that the condi- 
tion (45) is equivalent to beating the SQL bound. Then, 



as described in the previous section, entanglement is al- 
ways necessary and, in the symmetric case, also sufficient 
to violate the generalized uncertainty relation on the con- 
jugate parameters by means of the proposed setup. 



VII. REMARKS 

The estimation of the two conjugate parameters of 
a displacement operation is important both for appli- 
cations and fundamental reasons. Displacement opera- 
tions arc indeed ubiquitous in most of the quantum pro- 
tocols for CV systems. On the other hand, as we re- 
marked earlier, this estimation follows from the uncer- 
tainty relations, and thus from the foundational proper- 
ties of quantum mechanics. In this Letter we have pre- 
sented a measurement scheme which estimates accurately 
the two real parameters characterizing the unitary oper- 
ation of displacement in phase space, by using Gaussian 
entangled probe states and homodyne detections. We 
have further derived the ultimate quantum bounds on 
the multiparameter estimation for single and two-mode 
input Gaussian states, showing that our setup is opti- 
mal for a large range of parameter values characterizing 
the probe states. We have discussed the role of entan- 
glement, showing that in our setup its presence is always 
necessary, and in symmetric cases also sufficient, to beat 
the standard quantum limit achievable by using coherent 
input states and heterodyne detection. Finally we have 
analyzed in detail the relationship between our results 
and the generalized Heisenberg uncertainty relation for 
conjugate parameters. 
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